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Abstract. The fundamental measure density functional theory for hard spheres 
is generalized to binary mixtures of arbitrary positive and moderate negative non- 
additivity between unlike components. In bulk the theory predicts fluid- fluid phase 
separation into phases with different chemical compositions. The location of the 
accompanying critical point agrees well with previous results from simulations over 
a broad range of non-additivities and both for symmetric and highly asymmetric 
size ratios. Results for partial pair correlation functions show good agreement with 
simulation data. 
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Density-functional theory (DFT) is a powerful approach to study equilibrium 
properties of inhomogeneous systems, including dense liquids and solids of single- and 
multi-component substances p. Its practical applicability depends on the quality of 
the approximation to the central object, the (Helmholtz) excess free energy functional 
arising from the interparticle interactions. The specific model of additive hard sphere 
mixtures constitutes the reference system par excellence to describe mixtures governed 
by steric repulsion, and Rosenfeld's fundamental-measure theory (FMT) j2l El IH IHj is 
arguably the best available approximation to tackle inhomogeneous situations. A rapidly 
increasing number of applications to interesting physical problems can be witnessed jH]- 

The more general non-additive hard sphere mixture is defined through pair 
potentials between particles of species i and j, given as Vij{r) = oo for r < aij 
and otherwise, where r is the center-center distance between the two particles, 
and (Tjj is the distance of minimal approach between particles of species i and j. In 
a binary mixture non-additivity is measured conventionally through the parameter 
A = 2cri2/(crii + (722) — 1- The physics of non-additive hard sphere mixtures is 
considerably richer than that of the additive case. In particular the case of A > 
is striking, as small values of A are known to be already sufficient to induce stable 
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Figure 1. Illustration of the relevant length scales. The hard core interaction distances 
aij between pairs of particles of species ij = 11, 12 and 22 are related to radii through 
CTii — 2Ri, (Ji2 = Ri + Ri2 + R2 and 0-22 = 2_R2, respectively. The spheres of radii Ri 
and i?2 represent the weight functions Wa"^ and w^^^^ , respectively, and can be viewed 

(12) 

as "true" particle shapes. The sphere of radius R12 represents the kernel K^^ being a 
mere construct to generate the correct hard core distance (T12 between species 1 and 2. 

fluid-fluid demixing into phases with different chemical compositions (for recent studies 
see e.g. references [3 |H1 El CHI)- 

The treatment of general non-additivity is elusive in the FMT framework. The 
author is aware of successful studies only in four special cases: First, for the Asakura- 
Oosawa-Vrij (AOV) model jTH [T2], where species 1 represents colloidal hard spheres 
and species 2 (with (T22 = 0) represents non-interacting polymer coils with radius of 
gyration equal to (T12 — (o"ii/2), an excess free energy functional was given IT^. Second, 
a free energy functional for the Widom-Rowlinson (WR) model, where an = 022 = 0, 
was obtained jT3]. Third, the depletion potential between two big spheres immersed in 
a sea of smaller spheres was obtained through "Roth's trick" of working on the level 
of the one-body direct correlation functional [W\ ITT] . In this case the functional 
for the additive case is sufficient to obtain results, but the approach is limited to small 
concentration of big spheres. Fourth, in Lafuente and Cuesta's FMT for lattice hard 
core models, due the an odd-even effect of the particle sizes (measured in units of lattice 
constants), non-additivity of the size of one lattice spacing arises [TH]. This effect, 
however, is specific to lattice models and vanishes in the continuum limit. 

The aim of the present letter is to generalize FMT for hard spheres to the case 
of general positive and moderate negative non-additivity and arbitrary size asymmetry. 
The proposed extended framework accommodates, in the respective limits, the Rosenfeld 
functional for additive hard sphere mixtures [2], the DFT for the extreme non-additive 
AOV case ^Hj, and the exact virial expansion up to second order in densities. The 
structure of the theory, however, goes qualitatively beyond that of either limit. 
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The excess (over ideal) Helmholtz free energy functional is expressed as 



3 



^exc[pi,P2] = ksT I d^d^ J2 ^S?^(|x-x'|)<l>«^ ({n«(x)},{42)(x')}) , (1) 

a,f3=0 

where Pi(r) is the one-body density distributions of species i = 1,2 dependent on 
position r, ksT is the thermal energy, ^af3 for a, /? = 0, 1, 2, 3 is the free energy density 
depending on the sets of weighted densities {riu^} for i = 1,2, and the kernels K^^\r) 
are a means to control the range of non-locality between unlike components and depend 
solely on distance r. The weighted densities are built in the usual way 2J through 
convolution of the respective bare density profile with appropriate weight functions: 

j drpi{Y)wy{\y.-r\,Ri), i = 1,2, (2) 

where i/ = 0,1,2,3 labels the type of weight function, and Ri = aii/2 is the particle 
radius of species i = 1,2. The (fully scalar) Kierlik-Rosinberg form 1201 
Wu{r, R) is used in the following, as this renders the determination of the K^^^ (r) more 
straightforward. The Wu{r, R) are 

Wo = -6"{R - r)/(87r) + 6'{R - r)/(27rr), 
w, = 5'iR-r)/{8n), 
W2 = 6{R - r), 
W3 = Q{R-r), 

where R = Ri, the prime denotes the derivative w.r.t. the argument, S{-) is the Dirac 
distribution, and ©(■) is the step function. Alternatively, in Fourier space the weight 
functions are Wa{k, R) = Att J drwa{r, R) sin{kr)r / k and given as 

wq = c+ {kRs/2), 
w^ = {kRc + s)/{2k), 

(4) 

W2 = ATxRs/k, 



W3 = 47r(s - kRc)/k, 



with the abbreviations s = sm{kR) and c = cos{kR). The kernels K^^\r) in can 
be viewed as a/^-components of a second-rank tensor 

/ W3 W2 Wi Wq \ 



W2 wl Wq W_i 
Wi wl w'ii W-2 

\ Wo w-1 w-2 w_3 y 



(5) 



where indexing is such that the top row contains KqI^\ . . . , k'^^\ etc, and f distinguishes 
different elements. All K^^'^\r) possess a range of R12 = (J12 — (o"ii + cr22)/2, i.e. vanish 
for values of r beyond that distance (see figure [T] for an illustration of the length scales). 
The dimension of K^^^^ is (length) and hence the dimension of w^ is (length)"^"^. 
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The elements of K^^^) are defined, with R = R12 > 0, through Q, and furthermore 

w\ =6'{R-r), 

wl =6"{R-r)/{87r), 

wU = 6^'\R-r)/{64n'), 

w_, = 5"{R - r)/(27rr) - 5'^^\R - r)/(87r), 

w^2 = S^^\R - r)/(167rV) - d^'^^R - r)/{QATT^), 

w_3 = -6^^\R - r)/{87T^r) + 6^''\R - r)/{6A7r^), 

with the derivatives S^'^\x) = (PS{x) /dx"' for 7 = 3,4,5. Again we also give the Fourier 
space representation [being together with ^ also valid for R = R12 < 0], which reads 

wl = 4:Tc{kRc + s)/k, 

wl =c-{kRs/2), 

wl^ = -{k^Rc + 3ks) /{16n), 

= {k'^Rc - ks)/2, 

= -Pi?s/(167r), 

w_3 = ik^Rc-3k^s)/{16n). 

In order to express the dependence of the free energy density, $0/3 in equation ((T)), on 
the weighted densities Q we introduce ansatz functions Aa\ for species i = 1,2 that 
possess the dimension of (length)""^ and the order 7 in density (i.e. contain 7 factors 
nr^). Explicit expressions for the non- vanishing terms are 

<=n^o\ 4S=4M\ A^l=[n^^y/{2Arr), (8) 

41 = 4\ 42 = (40'/(8^), 4? = 4^ 43 = 1- (9) 

The excess free energy density is then constructed as 

= E E (4" + "f ') , (10) 

7=0 7'=0 

where ^l^'Jirj) = d^(pod{f])/dTi^ is the 7th derivative of the zero- dimensional excess free 
energy as a function of the average occupation number rj [3], ^odiv) = {^~v) ln(l— ?7)+77, 
and (p^iv) = V^odiv) 7 = 0. The specific form (fTIH) ensures both that the terms in 
the sum in possesses the correct dimension of (length)"^ and that the prefactor of 
fod in is of the total order 7 in densities, as is common in FMT. This completes 
the prescription for the functional; a full account of all details, also for multi-component 
mixtures and for lower spatial dimensionality, will be given elsewhere. 

Here we discuss some of the properties of the theory. For small densities it is 
straightforward to show that the correct virial expansion up to second order in densities 
is obtained, Fgxc — I d^rd^r' fij{\r — r'|)pj(r)pj(r')/2, where the Mayer functions. 
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fij{r) = exp{—Vij{r)/{kBT)) — 1, are recovered through 

/i2 = -E-i'^*<^*-?' U = -i:-^!^*w^L^ = l,2, (11) 

where * denotes the convolution, g{r) * h{r) = J d^r'g{r')h{r — r'). In the limit of an 
additive mixture, A — >• and hence R12 — > 0, one finds that Kapi^x) Q ii (3 ^ ?> — a 
and Ka{a-S){.x) 6{x) otherwise. This leads to a cancellation of one spatial integration 
in (0) and yields the Rosenfeld functional for a binary additive hard sphere mixture P] 
with radii Ri and R2- In the AOV limit, R2 ^ 0, one finds that — if a 7^ 0, and 
n^Q^ — s> P2 otherwise. The integration over x' in together with the kernel i^Q/3(|x— x'|) 
and the fact that the density no^''(x') = P2(x') appears linearly in $0/3, see A''^^ in (jHl), 
plays the same role that building weighted densities for the polymer species in the AOV 
case does. The resulting functional is equal to that for the AOV model [12]. However, 
in the WR limit, in contrast to ^1], terms higher than on the second virial level vanish. 
For A = — 1 the two species decouple, and J-'exc[pi, P2] = ^exc[pi] +^exc[P2] which is not 
obeyed by the present approximation, limiting its applicability to small values of A < 0. 

We next turn to an investigation of bulk properties of the theory. To assess 
structure, direct correlation functions can be obtained via 



{2)/| _ /|\ _ (5 PJ-' 2XC 



(12) 

pi ,p2=const 

(2) 

which can be shown to feature Percus-Yevick (PY) like behavior: clj {r > 0-12) = 0. 
Inverting the Ornstein-Zernike (OZ) relations permits to calculate partial structure 
factors, Sij{k), and partial pair correlation functions, gij{r). We have carried out Monte 
Carlo (MC) computer simulations in the canonical ensemble with 1024 particles and 10^ 
MC moves per particle; histograms of all distances between particles yield benchmark 
results for gij{r). We have chosen an intermediate size ratio of (T22/<^ii = 0.5 and have 
considered various values of A from —0.3 to 0.5 and a range of statepoints characterized 
by packing fractions, rji = npiafjG for i = 1,2. For A = 0, the current DFT reproduces 

the solution of the PY integral equation, as the functional reduces to the Rosenfeld 

(2) 

case, which is known to yield the same c^y (r) as the PY approximation. Results for the 
representative case A = 0.2 at two different statepoints are shown in figure |2l The core 
condition, gij{r < aij) = 0, is only approximately fulfilled, but the overall agreement 
between results from theory and simulation is quite remarkable. 

In principle one could envisage that this approach permits to study the depletion 
potential, V^g^pi"* (r) , between particles of species 1 being generated by the immersion into 
a "sea" of particles 2 through v2,^i{r) = —ksTln gii{r) for pi 0, and p2 = const. 
However, for the (relevant) case of small size ratios (e.g. c'"22/crii ~ 0.1, see |E1 IK!) 
already in both limits of additive hard spheres and the AOV model the results are only 
of rather moderate accuracy, underestimating the strength of the depletion attraction 
|13j . similar to results from the PY approximation. However, results from the present 
theory obtained through the OZ route (not shown) cross over smoothly between the 




' 
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Figure 2. Partial pair correlation functions, gij{r), between species ij — 11, 12 and 
22 (as indicated), as a function of the scaled distance t/cth, as obtained from the 
present DFT using the OZ route (dashed lines) and from MC simulation (solid lines). 
Results for 512 (322) are shifted upwards by one (two) units for clarity. Parameters are 
o'22/o'ii — 0.5, A — 0.3, 7]2 = ?7i/8 and r/i = 0.05 (lower), 0.1 (upper). For comparison, 
the theoretical critical point is located at 771 = 0.118, 772 = 0.0321. 



additive hard sphere case and the AOV case, similar to the correct behavior [THj 
Hence one can conclude that the pair structure predicted by the current DFT is similar 
to that of the PY approximation. This is a remarkable property, and one can anticipate 
test-particle calculations to yield superior results. 

Evaluating (P) at constant density fields yields an analytic expression for the 
bulk excess free energy for fiuid states, Fexc = ^exc[pi= const, p2= const]. The total 
Helmholtz free energy is then F = Fexc + ^bTV J2i=i 2 Pi\^^(PiM) ~ 1]' where Aj is the 
(irrelevant) de Broglie wavelength of species i, and V is the system volume. Via Taylor 
expanding Fexc in both densities one can show that it features the exact second virial 
coefficients (consistent with the correct incorporation of fij{r) on the second virial level) 
and also the exact third virial coefficients (see e.g. 0) provided 2(Ti2 > max(a"ii, cr22)- 

The fluid-fluid demixing spinodal can be obtained from (numerical) solution of 
\d'^{F/V)/dpidpj\ = 0, and the location of the critical point can be determined from 
minimizing one of the chemical potentials, fii or /i2, along the spinodal. Such results 
are compared in figure 01 to those from simulations for cii = (122, performed in the semi- 
grand ensemble by Jagannathan and Yethiraj jTU] and by Gozdz P, the latter study 
including a finite size analysis, for a variety of non-additivities ranging from A = 0.1 — 1. 
For the highly asymmetric case of (T22 = results from Gibbs ensemble simulations 
were obtained by Dijkstra [7j. For both size ratios the strong decrease of the total 
critical packing fraction with increasing values of A, as well as the overall functional 



Letter to the Editor 



7 



0.4 
0.3 
V 0.2 
0.1 


0.2 0.4 0.6 0.8 1 1.2 

A 

Figure 3. The total packing fraction at the critical point, rj^^f, where rjtot — iji + V2, 
for a non-additive binary hard sphere mixture as a function of the non-additivity 
parameter A. Shown are results from the present DFT (lines) and from simulations 
(symbols) for the symmetric case, tT22/o'ii = 1, by Gozdz [^j (filled squares) and by 
Jagannathan and Yethiraj ^01 (open circles), as well as for the highly asymmetric case 
of (T22/o'ii — 0.1 by Dijkstra 7 (crosses). 

dependence are very well described by the theory. However, the precise value at given 
A is underestimated. This behavior is not uncommon for mean-field like theories and is 
also present in the AOV case. A benefit of working on the level of the density functional 
is that the structure is consistent with the free energy. In figure |3] partial structure 
factors are shown for a range of values of A evaluated at the fiuid-fiuid critical point 
obtained from the free energy, and indeed Sij{k — > 0) ^ ±oo. 

In conclusion, having demonstrated the good accuracy of the predictions of the 
current theory for bulk fiuid properties of the non-additive hard sphere mixture, we are 
confident that it is well suited to study interesting and relevant interfacial situations, like 
the structure and tension of interfaces between demixed phases, wetting at substrates 
fIT\ and more. Note that any colloidal mixtures interacting with soft repulsive forces, 
as e.g. present in charge-stabilized dispersions, can be mapped (e.g. by the Barker- 
Henderson procedure) onto an effective non-additive hard sphere system. Hence one can 
anticipate experimental consequences of the structure and phase separation predicted 
by the present theory. The treatment of freezing |Hj requires additional contributions to 
the free energy functional j3 Hj . 

H. Lowen, R. Evans, R. Blaak and K. Jagannathan are thanked for useful comments. 
Support by the SFB TR6 of the DFG is acknowledged. This work is part of the research 
program of FOM, that is financially supported by the NWO. 
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Figure 4. Partial structure factors, Sij{k) for ij = 11,12,22 (as indicated), as a 
function of kan at the fluid-fluid critical point for size ratio (T22/o'ii =0.1 and non- 
additivity A = 0.2,0.5,1. The results for A = 0.5 (1) are shifted upwards by 5 (10) 
units for clarity. 
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